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ELECTROLYTIC HYDROGEN PRODUCTION - AN ANALYSIS AND REVIEW 


I. INTRODUCTION 


The potential advantages of hydrogen energy storage systems for 
electrical power generation is currently being assessed and recog- 
nized by many independent utilities, institutions, and agencies. 
These energy storage systems consist of three major subsystems: 
production, storage/transmission, and utilization. The current 
hydrogen scenarious agree that thorough understanding of a high 
efficiency, low cost production subsystem is first priority for 
this energy storage system. The favored near-term approach is 
the electrolytic production of hydrogen. An initial attempt to 
understand and model the characteristics of electrolysis cells is 
the subject of this report. 


This report attempts to analayze the inefficiencies encountered 
in electrolyzer operation. Following a discussion of the thermo- 
dynamics of electrolysis cell operation, current and near future 
technology of commercial units is reviewed. Information from 
manufacturers is then combined with some recent research efforts 
to form the basis of a mathematical representation of cell in- 
efficiency. The model so constructed is tested and thus is used 
to predict cell performance characteristics under various conditions. 


n. thermodynamics 

In the electrolysis of water, a current is passed through 
an aqueous solution to produce hydrogen and oxygen: 

R 2 0 + elec, energy — H 2 + 1/2 0 2 Cl) 

The electrical energy is thereby converted into chemical energy 
as hydrogen with a change in enthalpy at 25 G and 1 atm of 

£,H - 68,320 cal./g-mole 

The first law of thermodynamics for a steady- flow system is 

Q - w s = AH (2) 

in which 

Q " heat added to the system 
- W s t= useful work done by the system 

For the electrolysis unit, the work term is the electrical energy 
input to the cell, and is given by 

W e E (3) 

s 

in which 

tt = 2, the number of electrons transferred, 

9^ *= The Faraday constant = 96,500 Ccu/<mbs/gm equiv. 

«= 23,074 cal/volt - gm equiv. 

E *3 the electric potential applied to the cell, volts 


2 


( 4 ) 


Combining equations (2) and (3) gives 

E = A H - Q 
ni? 

For the present, interest lies in the reversible case (i.e., cell 
operation at thermodynamic equilibrium) for which equation (4) 
becomes 

E « AH - Qrev (5) 

n T 

Since an electrolysis cell operates isothermally, from thermodynamics, 


So that 


Qrev - TSS 


E m /iH-TSS 

rev ' nr 


The numerator of equation (7) is simply the change in the Gibbs 
free energy for the reaction, and at 25° C, 1 atm 

£ G = AH-T2JS = 56,690 cal/gmole 

The reversible cell potential is then 

E = a G (8) 


and at 25° C, 1 atm 

E ~ 56,690 = 1 *229 volts 

rev 2(23,074) 

Equation (8) may be combined with the Gibbs-Helraholtz equation 

-Ag + Ah = - t A6 (9) 

dr 


to give 


-E + AH = - T d Erev 
rev — -=r- — — 

nT dT 
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In order to integrate this equation, the variation of AH with 
temperature is needed,, Since AH is also a function of pressure, 
this equation is first written in terms of a standard state defined 
as 'pure components at 1 atm pressure . Denoting the standard reversible 
cell potential as Eo and the standard enthalpy change as AH° gives 

Eo A ~ T **0 for p = 1 atm (11) 

nr dT 

The dependence of Eo on temperature tan be formulated thermo- 
dynamically based on the observed heat of reaction, AH 0 . For the 

electrolysis reaction, the observed heat of reaction expressed 

as a voltage is 


( AH °/ur) = 1.480V @ 25° C 

If the heat capacity of water is taken as 1 cal/gm - °C and that of 
the gases Ibj and O2 as 7 cal/g-mole q -C (the value for an ideal diatomic 
gas) then the value of A IT 0 /n r can be expressed as 

A E* /n^ = 1.449 + 1.625 x 10" 4 T (12) 

in which the temperature, T, is to be expressed in °K. 

Combining equations (11) and (12), integrating and using the value 
of Eo at 25° C (1.229 volts) gives a relationship between Eo and 
temperature 

Eo = 1,449 + (1.877 - 1.625 In T) x 10“ 4 T 
a 0 ain with T in °K and (p - 1 atm) 


( 13 ) 


For ah electrolysis cell using a KOH solution as' the electrolyte 
the reversible hell potential varies with pressure and KOH concentration 
according to the Nerjjst equation 

E rev “ E o + ( R -' E / n 1 r ) ln PH2 (P02) (14) 

a v*o 

where 

E rev ** revers *kle ce H potential, volts 
E 

o *=> standard reversible cell potential, volts, 
ipas given in equation (13) 

R pt gas constant 

1?H2 =s pressure of produced, atm 

p 

02 *= pressure of 0^ produced, atm 

Hh 2 o " activity of water in solution 

The activity and pressure contributions to the reversible cell 
potential in equation (14) may be separated and simplified to 

E — Eo -r 3RT litF - RT In ^ H„0 (15) 

Tirr ^ 2 
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It- 


since the oxygen and hydrogen are generated at the same pressure 
(the cell operating pressure, 3? atm). With R=1.987 cal/gmol£,— °K 
this becomes 

E - 1.449 + (.646 In E + 1.877 - 1.625 In I - 
rev 

0.431 In O. H 2 0) x Itf 4 T (16) 

It remains to express the activity of water, ^0, in terms of 

KOH concentration. Using *the data of Costa and Grimes (ref. 1) 

with x equal to the weight fraction KOH in solution, 

♦ 

E ** 1.449 + r.646 InP + 1.877 + 2.7 x 2 
rev 

- 1.625 In Tj x 10" 4 T (17) 

Equation (17) is the desired expression for the reversible cell 
potential in terms of pressure, temperature and electrolyte con- 
centration. An actual electrolysis cell must operate at a cell 
potential of at least this value. It is instructive to return at 
this point to equations (2) and (3) to consider actual (not 
reversible) operation. 

Combining equations (2) and (3) gives 

Q r= AH “ l7E (18) 

For any set of conditions (V, T, X), the enthalpy change is fixed 
since enthalpy is a state property. The cell potential. 




E, however, depends upon cell design parameters, approaching E rev 
for a well-designed system. The actual cell potential required 
is always greater than the reversible cell potential corresponding 
to the same conditions of pressure, temperature and electrolyte 
concentration. ‘ for reversible operation, Q is a positive quantity 
over the range of interest since the enthalpy change is greater 
than the Gibb's free energy change, A. , which by equation (8) 
is the same as flT E. As the operation becomes more irreversible; 

the cell potential increases thereby decreasing Q. Eventually, 

■ 

a point is reached where Q vanishes and no heat need be added 
to the system. This is called thermoneutral operation, and by 
equation (18), corresponds to a cell potential 

E, - (19) 

* nT 

in which E is the thermoneutral cell potential. Vox even greater 
t 

irreversibility of operation, Q becomes negative. That is, in 
reversible operation heat is required, but cell inefficiency leads 
eventually to operation requiring cooling. 

Electrolysis cell thermal efficiency is always calculated 
based on a comparison of the actual cell potential with the cell 
potential for thermoneutral operation: 



The thermoneutral cell potential, Et, is used in preference to the 
reversible cell potential, E rev , since the former represents the 
total energy requirement of the, cell (4H), whereas the latter 
corresponds to on v that portion of the energy requirement which 
must be supplied as electrical work and not as heat ( AG). This 
is evident from equations (8), (18) and (19). 

Since water enters the system pure and not in combination 
with the potassium hydroxide, the AH value is not dependent upon 
KOH concentration. The effect of pressure on AH for this operation 
is rather small and neglecting this leads to the conclusion that 

and are equal. Using this in combination with equations 

(12) and (19) leads to 

E t » 1,449 + 1.625 x 10" 4 T (21) 

Equations (20) and (21) may he used to calculate the efficiency of 
an electrolysis cell. Since E can be less than Et, the efficiency 
can he greater than 1007, but operation in this range has yet to 
be demonstrated. 
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III. REVIEW OF CURREHT AMD FUTURE TECHNOLOGY 
III A. INTRODUCTION 

An electrolysis cell for the production of hydrogen 
from water consists of an electrolyte, usually an aqueous basic 
solution, circulating between two electrodes. A membrane prevents 
the generated hydrogen and oxygen gases from mixing, and separates 
the cell into anolyte and catholyte chambers, while allowing ion 
transfer between the too. In a basic solution, the cathode reaction 
is 

2 H 2 0 + 2 e ^ H 2 + 2 OH " 

and oxygen is produced at the anode 

2 0H‘ ^ 1/2 0, + HjC + 2 e 

The above expressions are overall reactions with no attempt being 
made to describe the still unsettled question of intermedia tes» 

The anode and cathode reactions combine to give a production of 
one mole of hydrogen plus one-half mole of oxygen from one mole 
of water. A basic solution is used in preference to pure water 
so as to increase the conductivity of the solution. This could 
also be accomplished using an acidic solution but the corrosion 
problems would be more severe. 

The earliest electrolysis cells were unipolar tank-type cells 
wherein a single cell could contain several electrodes, each of a 
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fixed polarity* The anodes and cathodes of the cell are separated 
by one or more membranes so as to prevent mixing of the hydrogen 
and oxygon gases produced. The most common present-day cell is 
of the bipolar type. .In this arrangement, each electrode serves 
a dual role as .anode and cathode. That is, the plate that serves 
as the anode for one cell also constitutes the cell wall and is the 
cathode for the adjacent cell. The cells are joined together to 
form a pack and pressed together between end electrodes which 
are charged by a direct current source. The result is a design 
physically resembling a filter press and these units are commonly 
referred to as bipolar filter press electrolyzers. Both tank-type 
and filter press designs are described by Mantel 1 (ref. 2). 

Considerable research was performed in the 1960's on defining 
the operating characteristics of water electrolyzers. The Allison 
Division of General Motors- (ref. 3) conducted a detailed investiga- 
tion. of the effects of current density, pressure and. electrolyte 
velocity on cell performance using 6N KOH. electrolyte at 75° C in a 
bipolar filter-press design. They also conducted testing and evalua 
tion of materials for electrodes and membranes, as part of a study 
t?hich encompassed all aspects of the design of an energy depot 
electrolysis system. Weight and. mobility were established as 
important criteria in this design. In contract, the study performed 


by Allis-Chalmers (ref. 4) considered low capital investment as 
the primary goal. This latter study was more restrictive but 
resulted in technology later applied commercially by Teledyne 
Isotopes. The thermodynamics and economics of the' Allis-Chalmers 
investigation was presented in the classic work by Costa and 
Grimes (ref. 1). 

The design and construction of industrial water electrolyzers 
was reviewed by Chapman (ref. 5) in 1965. More recently, Stuart 
(ref. 6) surveyed current operating performance of commercial 
units. Gregory's treatise on the hydrogen economy (ref, 7) provides 
a good review of current and advanced concepts and a more thorough 
review by the same author appears as part of an Q13R project report 
(ref, 8). Another recent review of water electrolysis was motivated 
by the need for an energy storage device (ref. 9). The latter is a 
particularly good source for a discussion of the many ways in which 
efficiency is defined in describing electrolyzer performance. A 
survey of water electrolyzers for oxygen generation has also been 
performed (re£. 10) for their importance in spacecraft life-support 
systems. 

The present summary of currently available and future generation 
electrolyzers has been compiled based on material found in the afore- 
mentioned reviews together with information received directly from a 


few commercial manufacturers of large-scale electrolytic hydrogen 
equipment (refs. 11, 12 , 13, 14, 15) „ 

III B. - ELECTROLYSIS CELL OPERATING PARAMETERS 

The thermodynamic reversible cell potential of an 
electrolysis cell using KOH electrolyte is given in equation (17) 
as a function of pressure, temperature and KOH concentration. 

Actual electrolysis cells require a greater cell potential (over voltage) 
due to the losses which occur in the cell. These losses consist of 
polarization at the electrodes, concentration polarization in the 
electrolyte solution, and ohmic losses in the cell. The losses in 
the cell depend upon the above mentioned operating parameters to- 
gether with such cell design parameters as electrode spacing and 
electrolyte velocity. The cell design parameters will be considered 
in a later section of this report. 

Electrolyzers are usually operated in the range 70-85° C with 
operation at the upper end of this range requiring somewhat higher 
pressures in order to minimize vaporization of the electrolyte. 

Both from a thermodynamic and kinetic viewpoint, higher operating 
temperatures are beneficial in reducing cell voltage requirements, 
and at least one manufacturer (ref. 12) is actively engaged in 


developing new materials which will allow operation at temperatures 
above 120° C while maintaining cell integrity. 

The effect of pressure on cell performance has been investigated 

(refs* 2, 3, 4, 16) with somewhat varying results. The reversible 

thermodynamic cell potential increases slightly as pressure is 

increased, but due to the reduction in bubble volume, there is a 

compensating reduction in the over-voltage. This point will be 

discussed in more detail later for the case of solid electrodes. 

» 

The use of porous electrodes to some extent offsets the effect of 
hubbies on cell performance but there are structural problems 
associated with using them at high pressure. It is to be noted 
that Teledyne Isotopes uses porous electrodes successfully at 
70 psig. 

An important operating parameter which does not appear in 
equation (17) is the current density. Basically, the greater the 
current density, the more inefficient the cell. A typical plot of 
cell efficiency vs current density appears below: 




The first steep portion of this curve is due to activation polar- 
Ization at the electrodes; the steady climbing is a result of 
ohmic losses in the cell; the second steep section indicates the 
approach to the limiting current density, that is, the current 
density corresponding to the maximum mass transfer rate consistent 
with the given cell configuration and operating conditions. 


The magnitude of the activation polarization term is dependent 
upon electrode material. ’For large industrial requirements, noble 
metals are too expensive. Nickel or nickel-plated steel is commonly 
used for the anode due to its availability, Id'? cost, corrosion- 
resistance in KOH solution and low activation over- potential char- 
acteristics. Electrode performance is improved by increasing the 
effective surface area^ thereby reducing the effective current 
density while maintaining cell size and current. All commercial 
manufacturers make use of this phenomenon in their design. 


The predominant effect of electrolyte concentration on <—11 
potential occurs due to the ohmic resistance of the electrolyte 
solution. Therefore, the cells using KOH solution as the electrolyte 
generally use 25-307. by weight KOH, thereby maximizing the conductivity 
of the solution. Electrolyte concentration also affects the rate of 
corrosion of cell components as well as the vapor pressure of the 
solution. This latter fact is exploited in a static feedwater 
design (ref. 17). 


1 / 
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All of the operating parameters discussed above affect 
cost and in each case, a variation to increase the efficiency (and 
hence decrease the operating cost) of an electrolysis cell also 

increases the capital cost (and therefore the fixed cost). De- 
pending upon cell configuration and design as well as economic 
factors t there is then an optimum choice of the operating para- 
meters which minimizes the sum of operating and fixed costs. 

The manner in which’ the operating parameters affect cost 
has been fairly well established. Higher temperatures increase 
efficiency thereby decreasing operating cost but the increased 
corrosion rate results in an increased maintenance and replacement 
cost. As previously described, increasing the* cell pressure 
increases the efficiency but again at the expense of increasing 
the capital cost because of the thicker wall construction and 
safety apparatus required. 

In order to discuss the relationship between current density 
and cost, it must be made explicit as to whether one is concerned 
with cost on an annual basis or cost per unit quantity of hydrogen 
produced. Consider, first of all, cost on an annual basis. As 
current density is increased, operating costs also increase due 
to a combination of increased current and increased cell potential. 


Capital costs don't change and the resulting sum of fixed and 
operating costs increases monotonically with increasing current 
density. However, if cost per unit quantity of hydrogen produced 
is considered, an increase in current density still causes an 
increase in the operating costs due to reduced efficiency, but 
since the hydrogen generation rate, is proportional to the current 
(and hence to the current density) the fixed cost is inversely 
proportional to current density. There is some value of current 
density which minimizes the sum of these costs. 

• 

III C. ELECTROLYSER CORPORATION LTD. 

The Stuart Cell manufactured by the Electrolyser 
Corporation is a tank type unit with steel electrodes, the anode 
being nickel plated. It operates at about 160° F and a pressure 
just exceeding atmospheric (lO'VG) using a 28% KOH solution as the 
electrolyte. The cell voltage is 2.04 v and the power consumption 
128 kwhr DC per thousand cubic feet of hydrogen. The size of an 
individual cell varies from the 4000 amp unit which produces 63.6 
ft"* H^/hr to a 22,000 amp unit producing 349.8 ft^ l^/hr. Larger 
requirements are handled by connecting individual cells in series. 
Electrolyser also markets packaged hydrogen generators producing 
anywhere from 20 to 1000 cubic feet of hydrogen per hour. 

The tank type electrolyzer has several advantages over the 
bipolar filter press design. It uses fewer, less expensive parts 
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and since construction is modular, repair or replacement of a single 

b 

cell can be accomplished by simply removing the defective cell via 
electrical connections. In a filter press type unit, the whole unit 
must be disassembled if a single cell diaphragm requires replacement. 

The tank type cell of course also has its disadvantages. Due to 
construction, it operates at rather low current densities, and generally 
requires greater floor space than a bipolar filter press unit. Tank 
type cells are also restricted with regard to operating temperature 
because of their size. 


The cell potential of 2.04 volts at an operating temperature of 
70°C for the Stuart cell can be converted into a thermal efficiency 
rating with the aid of equations (20) and (21). First, using equation (21) 


E * 1.505 v 
t 


So that 

-?7 = 1.505 x ICO = 74Z 
1 2.04 


The cell potential may also be compared to the thermodynamic reversible 
cell potential given by equation (17); 

E » 1.196 v 
rev 

calculated using an electrolyte concentration of 2S7. K011. The variation 
of cell potential with current density for the Stuart cell is shown in Fig. 1. 
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1X1 D e TELE DYNE ISOTOPES 


The Electra Cell manufactured by Teledyne is a bipolar filter 
press type using porous nickel electrodes and 25% KOH electrolyte at 
180°F. Teledyne manufactures these in three sizes, the smallest being 
hydrogen generators producing between 10.6 and 21.2 cubic feet of 
hydrogen per hour. Their intermediate size systems produce between 
10.6 and 4240 cubic feet of hydrogen per hour and their electrolysis 

plants produce 250 pounds or more of hydrogen per day with a nominal 

* 

cell potential of 1.84 v and a current density of 400 amps per square 
foot (or 4300 amps per square meter). Since the current density can 
vary depending upon the particular use to be made of the electro- 
lysis plant, Teledyne sizes their units using a computerized 
routine. The variation of required cell potential with current density 
for the Teledyne Electra Cell hydrogen plant is shown in Figure 1 along 
with a similar curve for an advanced cell which Teledyne expects to be 
commercial within the next three years. 

The cell potential of 1.84 volts at 180°F corresponds to 
a thermal efficiency of 

^ = 1.507 x 100 » 82% 

1.84 

where the thermoneutral cell potential of 1.507 volts was calculated 
using equation (21). This may be compared with the thermodynamic 
reversible cell potential given by equation (17) with 257. KOH at 70 psig 

E = 1.223 volts 
rev 
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Since Tele dyne builds hydrogen plants to order for a 
specific application and optimizes these plants using a computer 
program, the current density of the plant will depend upon economic 
factors such as the cost of electricity and the utilization of the 
plant. As an example, the cost of a 25 ton per day hydrogen plant 
with a 1007= utilization rate and a power cost of 10 mills per kilo- 
watt-hour is optimized at $6 million with current technology and 
$4.5 million for Che advanced cell. More expensive electricity would 
mandate more efficient operation (i.e., lower current density) which 
would tend to increase these figures. Electricity at 20 mills per 
kilowatt hour would increase these cost figures about twenty percent. 


Ill E. LURGI APPARATE - TECHNIK, GMBH 


Lurgi manufactures a high pressure Zdansky-Lonza cell 
which operates at 95° C and 30 atm" pressure with an electrolyte con- 
centration of 25% K0H. With a current density of 1900 amps per square 
meter, the cell potential is 1.84 volts. The Zdansky-Lonza cell is of 
the bipolar filter press type and uses nickel-plated wire gauze electrodes 
The variation of cell potential with current density for the Lurgi cell 
is shown in Figure 1. Lurgi is also developing a cell with a nominal 
current rating of 4000 amps per square meter in an effort to reduce 
capital costs. Lurgi estimates a cost reduction from 14.2 million dollars 
to 10.7 million dollars for a 25 ton per day hydrogen plant using the 



higher current density. These costs include the electrolyzer, sub- 
sidiary units, and the electrical rectification equipment* 


The thermal efficiency of the Lurgi high-pressure cell is 

= 1.509 x 100 = 827* 

Q 1.84 

where once agijn the thermoneutral cell potential was calculated using 
equation (21) . The thermodynamic reversible cell potential given by 
equation (17) with 257* KOH at 30 atm and 95° C is 

a 

E = 1.252 volts 
rev 

III F. BAMAG VERFA3REKSTECHNIK GMBtf 


The electrolysis unit marketed by Bamag is a bipolar unit 
operating at 80° C with a current density of 2500 amps per square meter. 
The pressure is essentially atmospheric with an electrolyte concentration 
of 267, KOH. The cell potential for the Bamag cell at the nominal current 
density is 1.92 volts and the variation of this quantity with current 
density is shown in Fig. 1* Bamag estimates an electrolysis plant 
producing 7.6 tons per day of hydrogen to cost 2.67 million dollars. 

The thermal efficiency of the Bamag electrolyzer is, with 


a thermoneutral cell potential of 1.506 volts. 


V - 1*50 6 x 100 = 787. 
1 1.92 


and the thermodynamic reversible cell potential is 


E = 1.186 volts 
rev 


Bamag is currently conducting research with the aim of producing larger , 
cheaper modules operating at a higher temperature. 


Ill G. NORSK HYDRO VERKSTEDER A-S 


This manufacturer supplies a bipolar filter-press electrolyzer 
which uses a 25% KOH solution at 80° C and a pressure of about 15 
inches WG. With a nominal current density of 1500 amps per square 
meter, the cell potential is 1.87 volts. The variation of this 
quantity with current density is shown in Fig. 1. 


The thermal efficiency of the Norsk Hydro electrolyzer is 


- 


1.506 x 100 « 817. 


1.87 


and the thermodynamic reversible cell potential is 


E * 1.185 volts 
rev 


Norsk Hydro estimates the cost of a 95 ton per day hydrogen plant 
as $39 million. No research is currently being conducted to improve 
efficiency or reduce costs. 
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Ill H. LIFE SYSTEMS, INC. 


The static feedwater concept developed by life Systems 
under NASA sponsorship was designed as part of a spacecraft waste 
reclamation system. According to Life Systems (ref. 17) this 
concept has terrestrial applications. 

Basically, static feedwater electrolysis is carried 
out by transferring the water to the cell as a vapor produced as 
a result of a vapor pressure difference between the electrolyte 
solution and the feedwater. The electrolyte solution is supported 
in a matrix between catalyzed porous electrodes. The electrolyte 
is a 35% KOH solution, more concentrated than conventional cells 
so as to increase the vapor pressure difference which causes transfer 
of water from the feed compartment to the electrolyte. At present, 
the Life Systems cell is designed to withstand pressures up to 
600 psi and temperatures up to 220 a F at a current density of 
600 amps per square foot (6460 amps per square meter) at a cell 
potential of 1.84 volts with projections into the near future of 
a cell operating at 2000 psi, 300° F with a current density of 
1500 amps per square foot (16,200 amps per square meter) but an 
increased cell potential of 2.02 volts. 

V The thermal efficiency of the present technology Life 

Systems cell is 

= 1.510 x 100 = 62% 

1 1.84 

and the thermodynamic reversible cell potential is 

E = 1.261 volts 
rev 

The variation in cell potential with current density for the Life 
Systems cell is shown in Fig. 1. 
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HI I. GENERAL ELECTRIC 


In contrast to the previous cells, 'ach of which 
uses a potassium hydroxide solution as the electrolyte, the 
electrolysis cell developed by GE uses a solid plastic sheet of 
perf luorinated sulfonic acid polymer as the electrolyte, This 
is coated on one side with a thin layer of platinum black to form 
the cathode and a similar thin layer of proprietary alloy catalyst 
is used for the anode. This system has been described in some 
recent papers by Titterington and others (ref. 18 and 19), 

• * 

The advantages claimed for this system are, operation 
at high pressure (up to 3000 psi), and high current density (greater 
than 1000 amps per square foot), and the elimination cf corrosive 
electrolyte which could carryover into downstream equipment. At a 
current density of 1000 amps per square foot (10800 amps per square 
meter) and a temperature of 180° F, the GE cell has a required cell 
potential of 1.85 volts. The thermal efficiency is then 

*7 » 1.506 x 100 = 81% 

1.85 

and the thermodynamic reversible cell potential is 

E = 1.176 volts 
rev 

The variation of cell potential with current density for the GE solid 
polymer electrolyte cell at 18G°F is shown in Fig. 1. 

The capital cost of the GE solid polymer electrolyte 

cell is currently rather high due to the high cost of the solid 

polymer electrolyte itself as well as the previous metals used for 
the electrodes. GE is conducting research to lower these costs by 
developing a thinner electrolyte or an alternative solid to be used 
as the electrolyte, GE is also attempting to increase the cell 
efficiency through an increase in the operating temperature. The 
problem here is one of preventing cell degradation at the elevated 


tempera tures . By 1985, GE expects to increase cell operating 
temperature from the current 220° F up to 300® F. • 

Another factor to be considered in cost comparisons 
is the operating pressure. The GE electrolysis cell is capable of 
operating at pressures up to 3000 psi thereby eliminating the need 
for a compressor when hydrogen is required at pressures, up to this 
value. On the other hand, this reduction in capital cost is. not 
accompanied by a corresponding decrease; in operating costs. Unlike 
-all the other electrolysis cells, the GE unit shows no increase in 
efficiency as the pressure is increased. In fact, the cell voltage 
increases with rising pressure due to the power required for,; com- 
pressing the hydrogen. 

Ill J. SUMMARY AND COMPARISON OF ELECTROLYSIS CELLS 

The efficiency of the various electrolysis cells is 
' shown in Fig. 1 in the form of a graph of the dependence of cell 
potential on current density. Typical design conditions and capital 
costs are summarized in Table 1. The capital costs represent only 
approximate figures. Accurate comparisons can only be made by 
establishing design criteria and requesting bids from manufacturers. 
Such items as degree of utilization, efficiency desired, and 
hydrogen pressure required will affect these costs. 


The efficiency values shown in Fig. 1 and in Table 1 
also deserve comment. As pointed out by Lurgi (ref. 11), the cell 
potential exhibits some fluctuation (^ 3%) and, according to Norsk 
Hydro (ref. 13), there is a net increase in this quantity with time 
of about 1 or 2 % per year. Some of the manufacturers quote cell 
potential values after one year of operation, but this is not true 
for all and consequently, there is some degree of uncertainty pre- 
sent when maxing cell voltage comparisons. 


IV MODELING OF ELECTROLYSIS CELL EFFICIENCY 


The difference between the actual cell voltage and the 
reversible cell potential is due to three factors: 

a) chemical polarization at the electrodes, 

b) concentration polarization in the electrolyte; and 

c) ohmic losses in the electrolyte and across the 
membrane. 

Each of these terms will be explored in some depth. 

IV A. CHEMICAL POLARIZATION 


Chemical polarization results from charge transfer inhibi- 
tion at the electrode and is determined by the catalytic activity of 
the electrodes and the surface roughness. According to Eisenberg 
(ref. 4), chemical polarization can be represented as 

AE chem = 01T/ «. nr ) A* < ) (22) 

where e< = transfer coefficient 

(jS 0 = exchange current density 
and st current density 

The exchange current density is 7.*elated to the height of the activation 
energy barrier, £ F* 

c5 = kT exp( * A F*/RT) (23) 

* 0 

where k is a constant which incorporates Avogadro's number, the Boltzma 
and Planck constants, the Faraday equivalent and the electrolyte concen 
tration. 


Equation (22) can be written in the more customary Tafel farm 

* 

AE chem “ 3 + b 1" $ W> 

in which the Tafel constants ”a n and u b n depend on temperature. 

According to equation (24 ) 3 the variation in cell potential with 
temperature, d ( & / dT, is dependent ;ipon electrolyte concen- 

tration and current density. 

An electrolysis cell gives rise to two chemical polar- 
ization terms, one each at the anode and cathode. These are referred 
to as the oxygen and hydrogen overvoltage respectively: 

^ E 02 = %2 + b 02 ln ^ 

^ E H2 = ^2 + b H2 ln ^ ^ 25 ' 

The values of the Tafel constants appearing in equation (25) have been 
reported by a number of investigators (ref. 20 to 23) with somewhat in- 
conclusive results. The Tafel slope for hydrogen, b^j is fairly well 
known at 25° C, and the constants for oxygen are greater than those for 
hydrogen. The influence of electrode material, electrolyte concentration 
and temperature on the Tafel constants is not very well known at this time 
although it is reporued that commercial manufacturers have considerably 
more proprietary information than what has been reported in the liter- 
ature. 

Based on information available at the present time, the following values 
have been chosen as representative 

£E q2 = .6 + .055 In $ (26) 

= - 3 + *°45 In $ 
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at 25°C for modern electrolysis cells with $ in amperes per square 
centimeter. The effect of temperature on cell performance has been 
reported (ref, 23) as follows e 

0^ - Nj « d < A E)/dT = -.0033 V/°C 

H 2 - Fe = d ( A E)dT = - .0025 V/°C 

at room temperature over a range of current density up to 2000 amperes 
per square meter. 

Combining this data with equations (22), (23), (24) and (26), gives a 
total chemical polarization loss as 

AE chem = S.iy + Q .{& 1 £ — O.SMt'T' (22) 

IV B, OHMIC LOSSES 


The Ohmic loss in an electrolysis cell is the sum of the 
iR terms due to anolyte, membrane and catholyte resistances. The 
conductivity of K0H and NaOH solutions exhibit maxima when plotted 
versus concentration at constant temperature, and consequently, the 
electrolyte concentration is often chosen to be at this value so as 
to minimize the ohmic loss. For KOH solutions at 25°C, the maximum 
conductivity occurs at 207= KOH, increasing to 357. KOH at 80° C. In 
cell designs where the gases produced bubble up through (or with) 
the electrolyte, the bubbles contribute to the ohmic loss. In 
addition, vapor bubbles (H^O) will also be produced further increas- 
ing the ohmic resistance. The volume rate of bubbles produced will 
depend upon pressure, temperature, electrolyte concentration and 
current density, and these factors together with electrolyte conduc- 
tivity, membrane resistance and cell design considerations determine 
the ohmic loss in the electrolysis cell. 
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Xn order to determine the contributions of the above . 
factors to the ohmic loss of an electrolysis cell., consider the model 
used by Funk and Thorpe (ref. 24): 

The cell consists of electrodes of width W and height H 
separated by chambers for ariolyte and catholyte flow with an internal 
membrane. Electrolyte enters the cell with an initial velocity in 
the cathode chamber and velocity V^q in the anode chamber. By the 
stoichiometry of the electrolysis reaction, twice as much hydrogen 
as oxygen is produced on a molal basis. Assuming ideal gas behavior, 
this 2:1 ratio is true on a volume basis as well. Xf the bubble 
velocities for 0^ and are the same, and if placing the 

membrane such that Jl\ = 2^2 then gives the same gas-liquid ratio 
in both chambers. The bubble volume, however, also includes the 
volume of water vapor generated, but since the partial pressure of 
H^O and the total gas pressure is the same on both sides of the 
membrane (assuming the gases are saturated with water vapor at the 
cell operating conditions), this 2:1 ratio is unaffected. Under the 
above assumptions, the resistivity of the electrolyte - gas mixture 
will be the same in both chambers* 

The ohmic loss in the cell may now be written as 

AE ohmic = ^ Rc (28) 
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where the effective cell resistance Rc is given by 

2 

Rc = 21 r _ X i$ (4i) + Rsa 

i,=l 1 


(29) 


where R s= membrane resistance, ohm- cm 
m 7 

= membrane resistivity (ohm-cm) * thickness (cm) 

^ electrolyte resistivity (ohm-cm) 

©£ = volume fraction of bubbles in the cell 

s= multiplier to obtain gas- liquid mixture resistance 
Subscripts: 1 for cathdlyte chamber, 2 for anolyte chamber 

Equation (29) contains, in addition to the aforementioned assumptions, 
the additional stipulation that the bubble boundary layer is formed 
quickly and fills up the entire chamber. Funk and Thorpe found this 
to be the case in their experimental investigation. Tobias (ref. 25) 
suggests calculating the resistance multiplier £( ok ) as 

£(©<)= 1 (30) 

1.5 

(1 -<*> 


whereas Ma shove ts 1 empirical equation (ref. 26) fits data in the range 
0 < < .74: 


f ( ) « 1 

1 - 1.78*+ ck 2 


(31) 


The value of ck will obviously vary in the direction of electrolyte 
flow, being zero at the cell inlet and increasing monotonically until 
the cell exit is reached. The value of at the cell exit may be 


calculated by noting that the rate at which gas leaves the cell is 
equal to the rate at which it is generated* Consequently, 


<*e = A$ 
nF 



p 


) 


where <Y is the volume fraction of bubbles in the electrolyte leaving 
e 

the cell, and O' is the slip ratio, defined as (gas velocity/liquid 
velocity). The value c< must be corrected in order to include the 
vapor in the bubbles. 


Assuming that the gas is saturated with vapor and letting p^s represent 
the vapor pressure of the solution at the cell operating temperature, 
the modified equation for the volume fraction of bubbles in the exiting 
electrolyte becomes 


d q ~ 


4 J T 


(32) 


A$RT + nf(P~^)<s-V,W-i 

Equation (32) holds for both the anolyte and catholyte chambers, and 
in fact, with the previously stated assumptions, yields the same 
value for each if the slip ratio is the same in each chamber. This 
is easily seen since Y\ = 2 for the cathode chamber and H = 4 for the 
anolyte chamber and - 2 $2, all other values being the 

same in the two chambers. In terms of volumetric electrolyte flaw 
rate per unit area of electrode surface^ V = W Vg/A, equation (32) 
becomes 





(33) 


The value of the void fraction to be placed into equation (31) must 
be some average over the electrode surface. For this one-dimensional 
problem, a simple arithmetic average gives 


^ = (1/2) 0< e 


(34) 


<*# 
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A better average is not easily obtained. The problem is that the 
current density varies in the flow direction, not the voltage* The 
bubble generation rate is, therefore, not uniform over the surface 
and the makeup of the voltage across the cell will vary with position. 

At the bottom (electrolyte inlet to the cell) the bubble resistance 
is low. As bubbles are formed, the ohmic loss increases and the 
current density decreases, thereby also affecting the chemical polar- 
ization contribution. 

The resistivity of the electrolyte solution, £ , is 

obtained from equivalent conductivity ( ) data as 

r f c 'imp 

in which N is the normality of the solution. The quantity is 

a function of N as well as temperature. Data from Perry (ref. 27) at 
18° C are corrected for temperature variations using 

* £ f 1S°C (36) 

1 + m (t - 18) 

. o„ 

where t = temperature, C 

m « temperature coefficient, .02 - .025 for bases. 

» 

Kirk and Othmer (ref. 28) suggest a value of .022 for m except for strong 
bases (such as KOH), but this value correlates well the data for 28% KQH 
(ref. 29). 

The final term in equation (29), the ohmic resistance of the membrane, 
depends on the membrane construction. Modern electrolysis cells use 
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an asbestos membrane which may be similar to the fuel cell asbestos 
studied by TRW (ref. 30) for which they report values as follows: 

2 

Rm = .24 ohm-cm for 20 mil thickness 

2 

= .36 ohm-cm for 60 mil thickness 

2 

D. Soltis (of Lewis) reports a value of .20 ohm-cm fox 10 mil asbestos in 
14NK0H. 


In summary, the ohmic loss in a water electrolysis cell may 
be calculated if the following quantities are known: the electrolyte 

concentration, pressure, temperature, electrolyte flow rate per unit 
electrode area, cell width, gas bubble slip ratio in the anolyte and 
catholyte chambers, current density and membrane resistance. The 
temperature and electrolyte concentration together determine the vapor 
pressure of the solution which has been correlated (ref. 31) for temp- 
eratures less than 25° C as 


P o 
*s 

» 1 - .055N (37) 

P ° 
r w 


where p o = vapor pressure of the electrolyte solution, and 
s 


p o 
r w 


vapor pressure of water at the same temperature. 


Data at higher temperatures (ref. 32) are also well correlated by this 
expression. Vapor pressure data for water may be obtained from standard 
sources. 


IV C CONCENTRATION POLARIZATION 


Concentration polarization results from concentration gradients 


32 


which exist in the neighborhood of the electrode and is small lor a 
cell with circulating electrolyte. This contribution to the cell 
voltage is usually either ignored or lumped with the chemical polar- 
ization terms. Since this term arises due to a resistance to mass 
transfer, anything that tends to reduce mass transfer resistance 
will reduce the effect of this term. Thus, concentration polarization 
may be reduced in any of the following ways: 


1) increasing the temperature 

2) increasing the electrolyte flow rate 

3) increasing the concentration of the species being 
transported; *i.e., use greater strength KOH solutions 

4) decreasing the spacing between electrodes 

5) decreasing the current density 

No data have been found in this area. Eisenberg (ref. 22) states that 


A E = RT 
cone —yr: 
n * 


In 




(38) 


where ^ is the limiting current density, a complex function of solution 
properties and hydrodynamic factors. For fully developed Poiseuille flow, 
Linton and Sherwood (ref. 33) give 


%■ af ( ¥■) 


'/l 


( t 


Ctrl 


fa 


where 


3 

C = ionic concentration, g moles/cm 
o 

V 0 = electrolyte velocity, cm/sec 

yf = thickness of electrolyte, cm 

2 

D = diffusion coefficient, cm /sec 
H - electrode length } cm. 
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The diffusion coefficient is given by Ibl (ref. 34) as 


D = KTVl * (40) 



where is the equivalent conductivity of the solution. 

The temperature dependence of — A_ can be obtained by combining equa- 
tions (35) and (36). Equations (38) through (40) may now be combined 
to estimate the concentration polarization. Note that /H in 

equation (39) is the electrolyte flow rate per unit electrode area, 

V in equation (33). 

The above calculations do not account for the presence of gas 
bubbles which hinder the diffusion of chemical species. This hindrance 
may be estimated by assuming that the active electrode area is reduced 
in proportion to the bubble fraction. The net result of this argument 
is to introduce a factor (1 ~0() into the right band portion of eq, (39). 


XV D CELL VOLTAGE SUMMARY 


The actual cell voltage may now be evaluated as 


E s= E + AE . + AE + AE . . 

rev chem cone ohmic 


(41) 


The reversible cell potential, E^ ev , is given by eq. (17); the chemical 
polarization term, i s shown in eq. (27). Ohmic loss is 

obtained by a complex analysis leading eventually to eq. (28); con- 
centration polarization is given by eq. (38) . 


The effects of operating variables and cell design on the 
actual cell voltage may now be stated in a qualitative way. Operating 
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temperature affects all the terms in eq. (41)* An increase in temp- 
erature decreases all but the last term in eq* (41). The ohmic loss , . 
depends on temperature through the void fraction and the resis- 

tivity r^. The negative effect of increasing void fraction due to 
increased temperature can be offset by operating at higher pressure 
and, possibly, with a higher strength KOH electrolyte solution to 
minimize vaporization, Xu any event, it would appear that operation 
at the maximum allowable temperature is desirable. Of course, there 
are also energy losses associated with cooling the gaseous products 

but the magnitude of these losses may be minimized by suitable heat 

% 

recovery systems. 

* 

As mentioned above, increasing the pressure has at least one 
positive impact for electrolysis cells; that is, it reduces the volume 
of the gas bubbles produced and also reduces vaporization of the solution 
which adds to the total gas flow for the system. There is, however, a 
beneficial aspect of vaporizing the water in solution. Namely, this 
vaporization process tends to maintain the cell operating temperature 
without the use of high electrolyte circulation rates which can result 
in the need for an appreciable amount of pumping power. In addition 
to the above, increasing the pressure increases the reversible cell 
potential and has an unknown effect on the chemical polarization and 
concentration polarization terms. 

The main effect of changing the concentration of the KOH 
electrolyte is on the ohmic loss in the cell. The resistivity of the 
solution passes through a minimum as KOH concentration is increased 
thereby suggesting a concentration level for cell operation, although 
the effect of KOH concentration on the concentration polarization has 
prompted consideration of the use of somewhat higher concentrations 
than those indicated by conductivity arguments (ref. 35) • In addition 


to these effects, there is the aforementioned effect of KOH concentration 
on vaporization of water from the electrolyte solution and the increase 
in E associated with increasing KOF concentration as ’dictated by 
thermodynamics and the Nernst equation. 

An increase in current density will not affect E but will 

rev 

increase each of the loss terms shown in eq. (41). The tradeoff be- 
tween operating efficiency and capital cost that results from this 
has already been discussed. 

Electrolyte velocity plays a role in the determination of 

the ohmic resistance of the cell and also the concentration polarization. 

* 

Increasing the velocity reduces the void fraction due to bubbles, helps 
to maintain a more uniform temperature profile and increases the limiting 
current of the cell; consequently, both ^ ^conc and ^ E ohmic are 
reduced. Selection of the optimum velocity is then a tradeoff between 
the above cell efficiency considerations and the power required for 
electrolyte circulation. 

Electrode design is a very important factor in the performance 
of a water electrolysis cell. Such items as construction material, 
surface preparation, electrode spacing and cell dimensions play a major 
role in the determination of polarization and ohmic losses in the cell. 
Electrode material significantly influences the Tafel constants for 
chemical polarization, particularly the M a M values. The choice of 
nickel and nickel-plated steel for the anode is made based on con- 
sideration of these Tafel constants vs. availability and cost. Erotic 
materials such as platinum can reduce these losses somewhat but not 
economically for large scale usage. For special situations, other 
factors could make the use of these materials desirable. 


Treating the electrode surface so as to provide a greater 
effective area also reduces chemical polarization losses. This can 
be done by using porous electrodes (ref. 11), ’’dimpling" or activating 
the electrode surface (ref. 13) or by employing finned electrodes 
(ref. 16) or electrode strips (ref. 36). In any event, the increased area 
acts to reduce the effective current density at the surface, thereby re- 
ducing the chemical polarization loss as given by the Tafel equation. 

The effects of electrode spacing and cell size come directly 
from the equations presented in the previous sections. Equations (38) 
and (39) indicate that concentration polarization decreases as the 
electrode size increases and the spacing between electrodes decreases. 
Starting with widely spaced electrodes, decreasing the spacing at 
first reduces ohmic losses due to the shorter distance as seen in 
equation (29). At some point, however, as the electrodes becoma 
closer, the effect of the gas bubbles becomes important and eventually 
dominates, causing an increase in the ohmic loss. These effects are 
shown in equations (30) and (32). There is then an optimum electrode 
Spacing which will result in a minimum ceil voltage. This optimum 
spacing may turn out to be less than can be effected by modern con- 
struction techniques (ref. 1). 

As indicated in the previous sections, rapid recirculation 
of the electrolyte will reduce both ohmic and concentration polarization 
losses. The ohmic resistance is decreased due to the sweeping away of 
the gas bubbles. This same sweeping effect also reduces the concentration 
gradients formed. The savings in cell voltage must be balanced against 
the increased pumping power requirement in an overall efficiency optimi- 
zation. 


Some of the methods for improving cell efficiency discussed 
above will result in increased capital cost requirements necessitating 
the establishment of tradeoffs between the two. One obvious area in 

which this occurs is in the selection of an operating pressure. A 

* . . 

higher pressure generally results in increased cell efficiency, but 

requires heavier and, hence, more costly construction, not only of the 

electrolysis module, but for much of the auxiliary equipment as well. 

There is also the safety aspect of handling hot lye under pressure 

which may dictate the use of more costly control systems, for these 

reasons, an economic study should consider fhe entire electrolysis 

plant. Similar arguments can lead eventually to the adoption of a 

systems study program which integrates the -electrolysis plant into 

* 

the environment of its use. Minimum cost of hydrogen produced may 
not be desirable, for example, when the hydrogen is being used as an 
energy storage medium for a large electrical power generating, facility. 
In such cases, electrolysis cell efficiency may be an overriding factor 
in producing electricity at the lowest possible price. Such considera- 
tions, however, are beyond the scope of this report. 

In summary, equations in che preceding sections allow the 
calculation of electrolysis cell voltage, which is inversely propor- 
tional to the cell efficiency, when the following quantities are known: 

Operating temperature 
Cell pressure 
Current density 

Electrolyte (KOH) concentration 

Tafel constants, oxygen and hydrogen sides 

Electrode- to-membrane spacings 

Membrane resistance 

Electrode size (height or diameter) 

Electrolyte flow rate per unit area of electrode 
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From the above list, several quantities may be optimized, whereas 
selection of some others is governed by considerations not related 
to cell efficiency.. The choice of an operating temperature, for 
example, cannot be made based on cell voltage since the latter is 
a monotonic function of temperature. Operation at the highest 
temperature which still allows the maintenance of cell integrity 
is then ..ndicated, There are, of course, tradeoffs in all consid- 
erations, and tell efficiency is just one consideration in electrolysis 

cell designo 


IV. E. Model Verification 

The algorithm presented in the preceding sections has been 
assembled as a computer simulation model for electrolysis cell poten- 
tial. This model was tested against both commercial and research 
cell data. Figure 3 shows the variation in cell potential with current 
density for three manufacturers, along with predictions obtained using 
the model. Although there are differences in slopes and intercepts 
between the commercial data and the model predictions, the model does 
show clearly the superiority of the Teledyne design. 

The model also correctly indicates the effects of pressure and temps 
temperature on eLectrolvsis cell potential as reported by investigators 
at Oklahoma State University (ref. 37), The values shown in Fige-es 
4 and 5 are in rather good agreement with the OSU data except at low 
temperatures. 


The differences between the data and the values predicted using 
the computer simulation model shown in Figures 3, 4 and 5 is most 
probably largely due to the lack of information concerning electrolyte 
flow rate* The cell potential is rather sensitive to variations in 
this parameter whose value is not reported by either commercial manu- 
facturers or academic researchers. In order to apply the simulation 
model, a simple equation was developed based on some of the cell 
potential data reported by investigators at Oklahoma State University; 

V *= .07 (4-.003T) (.003 + -) t 1,6 (42) 

P 

in which V is the electrolyte flow rate per unit electrode area, whose 
value is required in equation (33). Since the simulation model uses 
this empirical relation, it can -be expected to give the greatest 
agreement with data from which the equation was assembled. This is 
indeed the case. The best agreement occurs for the high temperature 
OSU data which was used to obtain equation (42). It is to be noted, 
however, that the agreement with commercial cell data shown in Figure 3 
is rather good when one considers that equation (42) was used in ob- 
taining the model predictions and this equation was developed inde- 
pendent of the commercial cell data. 


The simulation model was used to describe the effects of variations 
in pressure, temperature, current density and electrolyte concentration 
on the cell potential and its various components. For this purpose, 
a base case was selected and variations about this base case were 
examined. The base case operating parameters chosen were as follows: 


temperature ; 
pressure ; 
current density* 
electrolyte; 


375° K 
30 atm 

0.6 amps/cm^ 
6 N KOH 
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The effect of pressure on cell potential is shown in Figure 6. 

As the pressure is increased, the reversible cell potential increases 
due to the work of compression of the generated gases. On the other 
hand, the ohmic resistance drops due to a reduction in the volume of 
these gases. For the design shown in Figure 6, this results in a 
decrease in the electrolysis cell potential over the range shown, 
although in principle, an optimum pressure will be attained beyond which 
the increase in the reversible cell potential overshadows the re- 
duction in the ohmic loss. 

The corresponding picture for current density is shown in 
Figure 7. For lota current densities, the major cause of inefficiency 
is seen to be due to chemical polarization at the electrodes. As the 
current density is increased, the ohmic loss becomes increasingly 
important. 

The effect of temperature on cell performance is shown in 
Figure 8. Initially, increasing the temperature results in a decrease 
in cell potential due mainly to a reduction in chemical polarization. 
Eventually, a point is reached where the ohmic loss begins to increase. 
This is due to vaporization of the electrolyte solution itself as shown 
by the sharp increase noticed As the boiling point of the solution is 
approached. 

Figure 9 shows the effect of electrolyte concentration on the ohmic 
resistance of the solution. The minimum cell potential occurs almost 
exactly at the point of maximum conductivity, 

V. Electrolysis Cell Dynamics 

Electrolyzers have been suggested for the storage of energy derived' 
from either off-peak power, wind or solar energy through conversion into 
hydrogen. In such applications, the voltage supplied to the electrolyzer 
can be expected to vary over a rather wide range. Depending upon the 
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time scale of the voltage variations, the electrolyzer dynamics may 
be important, and in fact a favorable dynamic situation could be 
exploited through the use of an intentionally-varying voltage input. 

This has led to. suggestion of a pulsed mode of operation for elec~ 
trolyzers. 

It is not clear that intentional cyclic operation (e.g,, pulsed 
operation) will improve the operating characteristics of an electrolyzer. 
Alternate periods in which voltage is applied and removed may in fact 
have a detrimental effect on the system in that it may cause control 
problems for the gas removal system. The possibility of causing a 
decrease in gas purity also exists. For these reasons, Norsk (ref. 13) 
advises maintaining a base load equal to 257. of rated capacity at all 
times and going above this up to rated capacity or more when power is 
available. 

In pulsed electrolysis (i.e,, on-off operation), bubbles of gas 
produced during the "on" period have the opportunity to leave the 
system during the "off" portion of the cycle. The benefit achieved 
through the use of this pulsed electrical input lies then in the re- 
duced resistance, and consequently increased efficiency, during the 
"on" period as a result of the reduction in bubble volume. As a 
limiting case, the electrolysis cell model described earlier can be 
executed assuming zero bubble volume. This would give the best im- 
provement in performance that could be expected. The net reduction 
in cell potential is shown parametrically in Figure 10 for a particular 
cell design with operating parameters as listed. If the current density 
for the pulsed cell is the same as that for the normally operating cell, 
a reduction in cell potential of .2-. 3 volts can be achieved. This is a 
sizeable reduction but the comparison can be somewhat misleading. If 
the pulsed unit is the same size as the normally operating unit, then 
the pulsed cell will produce less hydrogen because it only operates 
during a part of the cycle. 
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The two inodes of operation may also be compared under conditions 
of equal hydrogen production rate. Since the pulsed electrolyzer only 
operates for a fraction of the time, it must operate at a higher 
current density during that time in order to compensate for the 
renaining time. Let B represent the fraction of the cycle for pulsed 
operation during which the electrolyzer is operating. Then if 0q is 
the current density for normal operation and 02 the current density 
for the pulsed mode, equal hydrogen production rates are obtained when 

B 02 ~ 01 

or (43) 

02 = 01 


The higher current density during pulsed operation will result in a 
decrease in cell efficiency. When this effect is combined with the 
improvement in efficiency as a result of removal of the gas bubbles, 
it is not clear whether or not pulsed operation is desirable. 

Suppose that the pulsed electrolyzer operates with equal time 
intervals for "on" and "off" sd that B - 1/2. Then equation (43) 
indicates that 02 = 2 0^. Defining voltage reduction once again as 
E pulsed - E normal, gives the voltage reduction curve labeled "equal 
hydrogen production rates" in Figure 10. In this case, the abscissa 
is 0q, the current density for the normally- opera ting cell and the 
value of 0i influences the choice between pulsed mode or normal operation. 
For the particular case shown, pulsed operation is desirable if 0^ is 
less than 0.5 amps/cm^ but is undesirable if 0^ is greater than 0.5. 

The dynamics of electrolysis cells should not be dismissed based 
on the preceding analysis. That is, the electrical energy source for 
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an electrolyzer may be generated by a varying power source such as 
photovoltaic or wind en< v.gy. Under such circumstances, the electrolyzer 
may be operated directly with this varying input, or it may be operated 
at a constant input through the use of a conventional power source when 
needed to supplement the primary energy input. Figure 10 indicates 
that this choice is controlled by electrolysis cell dynamics together 
with the dynamics of the energy source. 

VI Summary and Conclusions 

The analysis presented in this paper has led to a simple mathematical 
description of electrolyzer performance. This model is not claimed to 
be universally applicable and accurate - there are too many assumptions 
and unknown parameters for that to be true. However, the model does 
give a good representation of the effects of pressure, temperature and 
current density as is seen in Figures 3, 4 and 5, and this in turn 
endows the model with at least some credibility. As a next step, the 
model could be integrated into a multi parameter, optimization routine 
which would find the best designs and operating conditions subject to 
constraints such as temperature and pressure maxima. 

More importantly, the model indicates the relative magnitudes of 
each of the losses involved in electrolysis cell operation. For example, 
it was found that concentration polarization effects were negligible 
and for this reason, this term does not appear in Figures 6, 7, 8 and 9. 

In addition, the model helps to create a better understanding of these 
losses, and is used to predict the possible advantage of operating an 
electrolysis cell with a pulsed electrical input. 

Based on the present study, it appears that further research is 
needed to establish a more thorough understanding of the overvoltages 
in an electrolysis cell. This research should combine theoretical as- 
pects with experimental information. At the same time, research into 
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novel approaches to electrolyzer design and/or operation should be 

encouraged. Such research might include use of a pulsed electrical 

input or use of sunlight to reduce the voltage requirement of the cell 

through activation of the electrodes. Finally # a better understanding 

of cell dynamics is required in order to assess the feasibility of 

coupling an electrolysis cell to an unsteady power source. 

« 
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APPENDIX 


The Computer Simulation Model 

This program calculates the cell potential for an electrolysis cell of 
a given design with stated operating parameters. The design variables 
and operating parameters are the following: 

pressure = P, atm 

temperature = T, °K 

current density - PHI, amps/cm^ 

electrolyte (KOH) concentration = CONC, weight fraction, weight 

percent or normality 

slip ratio = SRA, SRC - ratio of bubble velocity to electrolyte 

velocity for anode and cathode chambers respectively 

flow rate ratio = HR, ratio of electrolyte flow rate in cathode 

chember to that in anode chamber 

cell width ratio = LR, ratio of width of cathode chamber to that 

of anode chamber 

electrolyte flow rate = W, cc's per second per square centimeter 

of electrode area 

cell width = L, cm 
cell height = H, cm 

2 

membrane resistance = RM, ohm-cm , = membrane resistivity (ohm-cm) 

*thickness (cm) 

Since in a given analysis some of these parameters may be unknown, the 
program uses Namelist input with a default option. In this type of in- 
put, the variables are assigned to a Namelist. For this program, the 
assignments are as follows: 

Namelist/OPVAR/PHI, P, T, CONC 
Namelist/RATIO/SRA, SRC, WR, LR 
Name list/CELL/W, L, H, RM 
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The data stack is composed of 3n+l cards in which n is the number 
of sets of data to be run as specified on the first data card using 
an 12 format. For each data set, there is one data card for each 
Namelist parameter list. These cards are typed with a M $" in column 
two followed immediately by the Namelist Name (OPVAR, CELL, or RATIO). 

The data is then listed in equation assignment form (e.g., p = 5, 

T = 400) with separating commas. Not all of the data need be specified. 
Those items not specified will remain at the values used in the previous 
data set. For the first data set, the default values are as follows: 

P = 1, T = 350, GONG = 6, PHI = 0.2 

SRA = SRC 1, WR = LR a 2 

L = 0.3, H = 100, KM = 0.1 

If W is not included in the specification of the first data set, it is 
calculated using equation (42) (V = W) . If it is desired to use this 
equation for subsequent data sets, a negative value must be assigned to W 
in the CELL Namelist. The program listing followed by a sample data set 
and the corresponding output completes this Appendix. 
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PROGRAM LISTING 


Main Program 


CCKMG*! /44/ 4<12> 
COMMON / 3 D/ DCEI 

:o«mov /: ;/ c i s i 

N 

DO ID I = l,N 
CALL T N T T 
C4LL PT%'TL 
CALL ° R N T R 
0CC10 CONTIVU^ 

GCC11 rOR‘MT H>1 
crcr 
END 


First Subroutine - Default Option 


C^OC* DATA 

C0 M MCV' /4A/ TtT »C°NC • FHT » RAT < 4 > » W r ;» » 

or 4L L 

DATA P.Tt CONC.PHl.lfAOIX »#: = 1*4I * 3IT. .e. .1*. 
1 2.*2.t-l.tC.3tlO"..r .10 / 

END 
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Next OUOIUULlUC - pa*-** **»*>*•»- 


SUG^OUTINE IN IT 

COMMON /AA/ P*T .CONCtPHTtSRA *SRC*WR»LR*WtL*H»nM 

COMMON / 33/ Vi A * W C * LI * L 2 • X » N 

REAL LR*L»Ll»L2»X*N 

NA MFL t ST /OPVAR/ P HI # P * T » C 0N0 

NAMELIS T /PATIO/ Z R 4 * 3 RC » WR * LR 

NAMEL T 5T /CELL/ W*L*H*RM 

WRITE (5.131) 

REA*' ( 5 t OP V A R ) 

WRITE (SiOPVARJ 

IF (CONC .LT. 0. .Of. CONC -5T. ICC.) 00 TO 130 
IF (CONC .ST. !. » 00 TO 111 
X = CONC 

00110 N = 1«> . *X*( X + l. ) 

GO TO 135 

DC 1 1 5 IF (CONC • GE • 15. J CC T" 12E 
N r CONC 

C01RP X r -0.5 + SCRT< C.2 r + (N/13. »> 
no TO 135 

001R5 X = 0 • G1 * CONC 
50 TO 11G 

0Q1 TO WRITE (&H3C) 

N - G. 

30 TO 12 U 

CC135 RE AC f 5 * R AT TO ) 

WRITE (S » RATIO! 

REAL I E * CELL > 

WRITE ( E * CELL ) 

W A r V/ ( WR+1 • I 
wc = W A * W R 
L 2 - L/tLR+l.J 
LI - L2*LR 

CG13G r ORMAT (44H ILLEGAL VALUE F0f CONCENTRATION* CN 
00131 r 0 R v A T (1H1) 

RETURN 
E NO 


ASSUMED ) 
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Subroutine to Calculate Cell Potential 


COROUTINE PTNTL 

COMMON /AAV P*TfCONC»PHI#SRA»SRC*WP*LR»W»L»H.RM 
COMMON / 3 3/ N 4 1 W C » Ll » L2 * X ♦ N 
COMMON /ZZt E * ER ; V * CCCHM f DL5HM C ^ I * OP CON 
PE AL L * L 1 » L 2 * L R * N » L * M * L A M IN 
C 

C T HIS SU3PCUTTNE EVALUATES THE CELL POTENTIAL 

c f OP '.HE DESIGN SPECIFIED 

C 

C THE REVERSIBLE CELL POTENTIAL - NEENST EQUATION 

00200 CREV r 1.443 ♦ IC.C4£»4L0G I PI +2. 7* X **2*1 .377-1 *S2E *ALOG( Tn*T*10 
1 * ( - a ) 
c 

C CHEMICAL POLARtZATlOV VIA THE TAFEL E3JATICN 

C MODIFIED F OP TEmo?R4TuRE VARIATION 

OECHM = 5.13 + C.!*AwOG|PHI> - 0.2*AL0 r m 
IF < D r CHM .LT. 0. ) CECHM=0. 


C VAPOR PRESSURE CCPR:LAT TONS 

PW0-1P.**! 7.117 - < ''222. /T ) J 
PSC = PWC*U.-0.t^5*M»/l4.7 
IF <PEu .LT. PI 30 TO 2 If 
WRITE (Ht 2341 PSO 
P r 1 . 2 * PSD 

00215 IF f VI .CT. P.| CO T" 220 

w ~ d • 07* 1 4 • -c . nc’ T * r i * i n. rc ?+ 1 1 ■ /Pi i * i Phi * *i . e i 
W4 - W/fWR+l.l 
WC r WA*rfR 

c 

C EU33LC VP^JML CAL^ULATICX 

00220 DC - 2 25 2 4 *5R C* WC * ( P -P EL ) / ( T* PHI ) 

CA r ?.*SS*SPA/IS!?C*WR) 

a*.pwc = l./tac+i.t 

ALPHN r 1 » / { <54 ♦ 1*1 


C 

C 


C 

w 

002 20 

C 

c 


c 


RESISTANCE MULTIPLIER C7RREL A ^ON 
FC = l./l 1 .-ALP RC!** 1.5 
~C = U.+2 .*f:>/3. 

FA r I ./ 1 1 4-ALPHN ) **! .5 
"fl - f 1.+2. *FA)/3. 


j^DUCIBILnY OPTOE 

PAfiB IS POOR 


CC PEEL ATI ON FO® EQUIVALENT CONDUCTIVITY 
LAM - 223. *1D .*♦ I-NV 15*! 

RESIST IV TTY OF KG E4 WITH TEMPERATURE CORRECTION 

PFA = 1000./ IN*L4M1 

RF = RF4/ C 1. + 0.022*1 T— 291 * 1 J 
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Cell Potential Calculation - Continued 


e OHMIC LOSS CALCULATIONS 

DEOHM* 4 > ' PHI *RK 
RE OHM ( 2 ) s L?*FA*RF*?MI 
DECWMf 31 r Ll»-C*RF*PHI 

3C0HM«l» = SE 0 HM ( ? 1 ♦' DEOHMdl ♦ DC OH Ml 4 1 

c : ; ' 

C E3UIVALENT CONDUCTIVITY 

C AND DIFFUSION COEFFICIENT CALCULATIONS 

00240 L AMIN r 1G0C./(N*RFI 

D = t, AED *10 • ** ( -ID ) *T *LANI N 
C 

C LIMITTNP CURRENT *'ENSITY CORRELATIONS 

n HlL T M - S?1.4*N*V4**t l./ 3.TMD/L2> **( 2,71. 1 

C CHECK ON ASSIGN ED CURRENT DENSITY 

IF (PHI .lT. PHI. IN ) C n ID 7 r C 
WRITE (C»2Sr» PHI*P HIL I M 
DEC CN s 0. 

GO T 0 2<>C 
C 

C CONCENTRATION POLARIZATION 

002E0 S r ?.1?*10.**( -S)*T 

D.CCN - 3* ALOG IPHTLTM/ C°HXLI m -PHI ) I 
C 

C CELL POTENTIAL 

00290 E = EREV ♦ OECHM ♦ DEGHMt! I ♦ DE CON 


00280 FORMAT (42H ***** CURRENT DENSITY "YCEEPS LIMIT ♦•*♦♦» 
1/ 5H P HT = *F1C.4 / 3H °HILIM r , ^ 1 0 • J 
002P4 FORM A T (79 HO T HE VAPOR pPESSJ^E OF THE ELECTROLYTE EYCEE 
IS MED OPERATING PRESSURE / 44H PRESSURE REASSIGNED AS 1 
2RE PSD = • FI 0.4* ATM ) 

return 

END 
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Subroutine to Print Out Results 


■3U.9.*miItN:£ PRNTR 

COMMON / A 4/ P*T*CONC*FHT»SR® • SRC • WP* LR * W * L #H * "*M 
COMMON / 33/ WA*WC*Ll *L2*X#N 
COMMON /CC/ F*lREV#P£CrtM*DE3MMl4T*?EC0N 
NEAL LRtL *L1*L 7* X*N 
C 

C PRINT OUT DESIGN CPI GIF T&4 7 IONS 

0030* WRITE (G.401) 

WRIT* (1*402) P*T*PHl#« 

C 

C *RTMT OUT CELL OUT A 

WRITE (5*433) 

MPT 1 ' G (S»4P4t W * !_ * M * 7 M 

A 

C SPECIFICS 

WRITE (5*405 I 
WRITE (G*40G> WC#WA 
W R I T E ! u * 4C7 ) HtL2 
WRIT*- *5*403) SRC*SRA 

C 

C RESULTS 

0C33CT WRITE (5*412) E 

WRI^E (5*413) EREV 

WRITE (5.414) DE:mm 

W R I T G (S*41S) ( OE*>rP' ( I ) * Is 1 *4 ) 

WRITE (E.415) OECRN 

QG4C1 F ORM A T I 42Hl***»***** DESIGN SPECIFIC A f TONS ********* ) 

00402 *OR M fiT (lunpRESSURT “ * FJD.3* 4H 4 Ty / 14 if TEVFERATURE r * F10*1 
1* GM REG K / 19H CURRENT {JCNCITY ^ * Fl r '-' > * AMF/CDC U / 

2214 KOH CONCENTRATtcv - # "IP. 2* TH NORMAL ) 

00403 FORMAT (2CHC**** "ELL DATA **♦* I 

00404 *OPMflT (PAhOElICTPClVTE F|.:W RATE t * FlO. 4* 12H CC/SEC-SDCM / 

11 3H CFlL WIDTH - » F1L.2* 7H CM / ?GH r LECT POPE HEIGHT z * ^IF.1* 

2 34 C v / 2EH MEva^A^E RESIST ANC r i » FlO.' 1 ’* 9H OHM-CMS 3 ) 

00405 r 0P M A T ( 2 0 HO * * * * SPECIFTCS **** / 23X* !7HC4 ThcDE A»'DDC ) 

0040G FORMAT (174 ELECT^OlYTE RATE • 7X* 2 r 1 C. 4 ) 

00407 r 0RMA7 (11H CELL WIDTH ♦ 9X* 2F1G.4I 

00408 r 0RM4 T (114 SLTP * 3X* 2F1C.4I 

00412 r DPMA T I 1H0 / 23HC********* RESULTS ********* / IHOtlOX* 

1214T4E CE^L POTENTIAL IS » r 10,4» *H VCLTS / 

223H0CPNSTRUCTED AG ALLOWS ) 

0041** POR«AT ( 1 HO • 1QX* IRREVERSIBLE CELL POTENTIAL = * r 10.4* 

1 5H VOlTS ) 

0C414 PORMAT ( 1 .40 » 10X* Z’HCHIMICAL P0LARI7A TTON r * FlO. 4* 6M VOLTS ) 
00415 FORMAT ( 1 HO * IPX# 2QH0HMIC LOSS IN CELL - • FIG. 4* GH VOLTS / 

1194 MADE UP AS FOLLOWS / 13H 4NC3G SIDE • F1C.4 / 

213H CATHODE SIDE *F!C.4 / 9H MEMBRANE # F1C.4) 

004 IE r 0R«4T (IHOtlOX. *>8 hCCNCENTR ATION POLARIZATION - , F1C.4*EH VOLTS) 
RETURN 
END 


REPEOBUCIBILffY OP THE 

ORIGINAL PAGB IS POOR 
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Sample Data Set 

01 

40PVAK 

$*477 0 $£VZ> 

4cbll 


Output for this Data 
/ 


********* DESIGN SPECIFICATIONS ********* 


PRESSURE - 1.00G ATM 

TEMPERATURE =. 350.0 DEE K 

CURRENT DENSITY r .20 4MP/S0CM 

KOH CONCENTRATION r ".00 NORMAL 


CELL DATA ♦♦** 

ELECTROLYTE FLOW RATE = .0159 CC/SEC-SSCM 

CELL WIDTH = .30 CM 

ELECTRODE HEIGHT = 1GC.0 CM 

MEMBRANE RESISTANCE = . 100 CHM-C‘4S3 


•*** SPECIFICS **♦* 

ELECTROLYTE RATE 
CELL WIDTH 
SLIP RATIO 


CATHODE ANODE 

.CUE* .C0S3 

.2000 .1000 

l.coro i.gopo 


•*****••* RESULTS *•••**•** 

THE CELL POTENTIAL IS 1.0140 VOLTS 

CONSTRUCTED AS FOLLOWS 

REVERSIBLE CELl POTENTIAL = 1.2381 VOLTS 

CHEMICAL POL ARITAT T ON = .3427 VOLTS 

OHMIC LOSS IN CELL - .78?G V^LTS 

HADE UP AS FOLLOWS 
ANODE SIDE .1212 

CATHODE SIDE .2424 

MEMBRANE .0200 


CONCENTRATION POL API? AT ION 


0003 VOLTS 


Table I 


Manufacturer 

Electrolyser 

Teledyne 

Lurgl 

Batnag 

Norsk Hydro 
Life System 


Summary of Electrolysis Cell Performance 


* 

atm 

X 

°c 

2 

amp/cm 

nun v' 
% 

1 

70 

.16 

28 

6.8 

82 

A3 

25 

30 

95 

.19 

25 

1 

80 

.25 

* 

* 

26 

1 

80 

.15 

25 

1 

93 

.65 

35 

200 

82 

1.00 



E Size Cost 


volts 

% 

tons /day 

10 6 § 

2.04 

74 

- 

- 

1.8/ 

82 

25 

6 

1.84 

82 

35.7 

/y,z 

1.92 

78 

7M 

2 47 

1.87 

81 

95 

3S.«? 

1.84 

82 

. 



G. E 


1.85 


81 


EFFECT OF CURRENT DENSITY ON ELECTROLYSIS CELL POTENTIAL 
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MODEL VERIFICATION 
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Figure 5 

MODEL VERIFICATION USING RESEARCH CELL DATA 
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Figure 9 


EFFECT OF ELECTROLYTE CONCENTRATION ON CELL POTENTIAL 
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